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We investigate Gaussian quantum states in view of their exceptional role within the space of all 
continuous variables states. A general method for deriving extremality results is provided and ap- 
plied to entanglement measures, secret key distillation and the classical capacity of Bosonic quantum 
channels. We prove that for every given covariance matrix the distillable secret key rate and the 
entanglement, if measured appropriately, are minimized by Gaussian states. This result leads to 
a clearer picture of the validity of frequently made Gaussian approximations. Moreover, it implies 
that Gaussian encodings are optimal for the transmission of classical information through Bosonic 
channels, if the capacity is additive. 



States with a Gaussian Wigner distribution, so called 
Gaussian states, appear naturally in every quantum 
system which can be described or approximated by a 
quadratic Bosonic Hamiltonian. They are ubiquitous in 
quantum optics as well as in the description of atomic en- 
sembles, ion traps or nano-mechanical oscillators. More- 
over, Gaussian states became the core of quantum infor- 
mation theory with continuous variables. 

Besides their practical relevance, Gaussian states play 
an exceptional role with respect to many of their theoreti- 
cal properties. A particular property of Gaussian states is 
that they tend to be extremal within all continuous vari- 
able states if one imposes constraints on the covariance 
matrix (CM). The best known example of that kind is the 
extremality with respect to the entropy: within all states 
having a given CM, Gaussian states attain the maximum 
von Neumann entropy (cf.Q). Similar extremality prop- 
erties have recently been shown for the mutual informa- 
tion and conditional entropies 0, . 

In this work we prove extremality results for Gaus- 
sian states with respect to entanglement measures, secret 
key rates and the classical capacity of Bosonic quantum 
channels. These findings are based on a general method, 
which exploits the central limit theorem as a active and 
local Gaussification operation. Our main focus lies on 
the entanglement, which will serve as a showcase for the 
general procedure. We prove that for any given CM the 
entanglement, if measured in an appropriate way, is lower 
bounded by that of a Gaussian state. The same result 
is shown to hold true for many other quantities like the 
distillable randomness and the secret key rate. This re- 
sult not only emphasizes the exceptional role of Gaus- 
sian states, it also leads to a clearer picture of the va- 
lidity of frequently made Gaussian approximations. In 
practice, states deviate from exact Gaussians and their 
precise nature remains mostly unknown. Nevertheless, 
the CM can typically be determined, e.g., by homodyne 
detection, and one is tempted to calculate the amount 
of entanglement, or other quantities, from the CM un- 
der the assumption that the state is Gaussian. The de- 
rived extremality of Gaussian states now justifies this 



approach as it excludes an overestimation of the desired 
quantity, even in cases where the actual state is highly 
non-Gaussian. In this sense one stays on the save side 
when a priori assuming the state to be Gaussian. We 
will see, however, that some care is in order, since the 
extremality property with respect to the entanglement 
turns out to depend on the chosen entanglement mea- 
sure. 

Before we derive the main results we will briefly re- 
call the basic notions. Consider a Bosonic system of 
N modes characterized by N pairs of canonical oper- 
ators (Qi,Pi, . . . ,Qn, Pn) ='■ R or equivalently by A 
Bosonic annihilation operators aj = (Qj + iPj)/y2, For 
any density operator p of the system we define a vec- 
tor of means with components dj = tr [pRj], a CM 
Tfc( = tr [p{Rk — dk, Ri — di} + ] and introduce a charac- 
teristic function = tr [pexp [££ • R]], £ € R 2N . The 
latter is the Fourier transform of the Wigner function 
and it thus completely characterizes the state For 
Gaussian states the characteristic function has the form 

x(0= e*<-i* r «, (1) 

such that they are entirely specified by d and T lead- 
ing to a complete description within a finite dimensional 
phase space M. 2N . The underlying Hilbert space H is, 
however, infinite dimensional and we will denote the set 
of all bounded linear operators on Ti by B{TL). Note that 
the following results also apply to finite dimensional sys- 
tems by simply embedding C d into Tt. 

Our results are based on a non-commutative central 
limit theorem, as discussed in 0, 0j an d operator- 
topology arguments from 0, Q . We will first state the 
main ingredient as a general Lemma and then discuss 
its applications to quantum information theory. Readers 
who are mainly interested in the applications may skip 
the proof of the Lemma. 

Lemma 1 Let f : B(H® N ) — ► R he a continuous func- 
tional, which is strongly super- additive and invariant un- 
der local unitames f(U® N pW® N ) = f(p). Then for 
every density operator p describing an A -partite system 



2 



with finite first and second moments, we have that 

f(p) > f(pa) , 



(2) 



where po is the Gaussian states with the same first and 
second moments as p. 

Let us first remark on the requirements in Lem.^ 
Continuity is understood in trace-norm, i.e., \\p^ — 
p\\i — > should imply f(p^ n ') — > f(p)- In fact, lower 
semi-continuity suffices, and we may restrict the domain 
of / to an appropriate compact subset of density opera- 
tors like those satisfying an energy constraint. The lat- 
ter typically restores continuity for functionals which are 
known to be continuous in the finite dimensional case. 
Strong super- additivity means that given a state p acting 
on (TCai ® 7~tA 2 ) ® 0~Lbi ® 7~Lb 2 ) with restrictions p t to 
Ti-Ai ® Kflj, then /(p) > /(pi) + /(/02), with equality if 
p = pi ® P2- The latter is referred to as additivity and 
both properties are analogously defined for more than 
two parties. 

Proof: The main idea of the proof is covered by the 
following equation: 
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In the first line we use additivity of / and set p = 
U® N p® n U^® N , where U is a suitably chosen local uni- 
tary which acts on n copies of the state. In the second 
line we first exploit strong super-additivity in order to 
bound /(p) from below by the sum over all reduced den- 
sity operators pV Then we argue by the central limit 
theorem and a special choice of £/, that each of these re- 
duced states pk converges to the Gaussian pq in the limit 
n — ► 00. 

This idea is now made rigorous in two steps. First 
we prove that each characteristic function \k converges 
pointwise to the corresponding Gaussian Xd an d then we 
argue that this implies trace-norm convergence on the 
level of density operators. To simplify matters we will 
w.l.o.g. assume that p has vanishing first moments, i.e., 
dj = 0. The general case is then obtained by applying 
local displacements, which by assumption will not change 
the value of /. 

Let us begin with specifying the local unitary U as 
a passive symplectic operation acting on the canonical 
operators on site a £ {1, ... , N} as 



la,k 
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—j= Pad , (5) 



(I 1 \® m 

with H = I j _i J being a Hadamard matrix and 
n — 2™. Physically U corresponds to an array of 50:50 



beam splitters and half-wave plates. Note that H has 
only entries ±1, so that we can partition the sum 



1 / "+ n - 

Qa,k = j — I ^ ^ QaA ^ ^ Qa, 
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(6) 



l:H kl = l 



j:H kj =-l 



and similarly for P a ^- Here n + = n — ?i_ is the number 
of ones in the k'th row of H. Note that either n + = n in 
the first row, or n + = n/2 in all other rows. 

The characteristic function Xk of the reduced density 
operator p^ is then given by 

/ N 



Xk{q,p) = tr 



QaQot.k ~t~ VolPql, k 
\ a = l , 



(7) 
(8) 



where x is the characteristic function of p and £ = (g, p) . 

Following |g| we introduce a function g : E — > C by 
9i x ) — x( x £)i wn i cn is a classical characteristic func- 
tion, i.e., the Fourier transform of a classical probability 
distribution with second moment £ T r£/2. To see this 
note that x is the Fourier transform of the Wigner func- 
tion and recall that every one-dimensional marginal of a 
Wigner function (in particular the one corresponding to 
the direction £) is an admissible probability distribution. 
Characteristic functions are continuous at the origin, sat- 
isfy 3(0) = 1, \g(x)\ < 1, and in the case of finite second 
moments we can expand up to second order |lCt , such 
that 

g( X ) = l~^R X l + (x 2 ). (9) 

Pointwise convergence Xk — y Xg follows then from Eq. JHJ) 
together with setting x = 1 in 
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For the remaining part we can combine the argumenta- 
tions in Refs. 0, E| ■ In it was proven that pointwise 
convergence of the characteristic functions implies con- 
vergence of the respective density operators within the 
weak operator topology. The latter was, however, shown 
to be equivalent to the trace norm topology on density 
operators in Ref.[9j. 

□ 

A simple application of Lem.^is the rederivation of 
the maximum entropy principle by setting / equal to 
minus the von Neumann entropy S(p) = — tr[plogp]. 
Similarly, in the bipartite case with N — Na + Nb and 
f(p) = S(pa) — S(p) we recover the recentlyproven ex- 
tremality result for the conditional entropy jjj for which 
strong super-additivity is an immediate consequence of 
the strong sub-additivity inequality for the entropy. 
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ENTANGLEMENT MEASURES 

Grouping the N tensor factors in Lem.^into M < N 
parties and exploiting that every entanglement measure 
is by definition invariant under local unitaries, yields the 
following. 

Proposition 1 Let E be a continuous entanglement 
measure which is strongly super-additive. Then for ev- 
ery density operator p describing an M -partite system 
with finite CM (and arbitrary, finite, number of modes 
per site), we have that any Gaussian state pa with the 
same CM provides a lower bound E(pc) < E(p). 

Most of the entanglement theory developed so far is 
devoted to bipartite systems. Whereas for pure states 
there is an essentially unique measure of entanglement, 
the von Neumann entropy of the reduced state, there are 
various different entanglement measures for mixed states 
|l2| . Two of them have a clear operational meaning: 
the distillable entanglement Ed quantifies the amount of 
pure state entanglement that can asymptotically be ex- 
tracted by means of local operations and classical com- 
munication (LOCC), and the entanglement of formation 
Ep (resp. its regularized form, the entanglement cost E c ) 
measures the pure state entanglement required in order 
to prepare the state. Among all other measures, the log- 
arithmic negativity Ejy is the most popular one, as it is 
comparatively easy to calculate Let us now discuss 
the consequences of Prop. for some entanglement mea- 
sures: 

(i) Distillable entanglement: Ed is additive (due to the 
asymptotic definition) and strongly super-additive (since 
restricted protocols lead to smaller rates). It was shown 
to be continuous in the interior of the set of distillable 
states 0|- As distillable Gaussian states are always in 
the interior [l^, Ed fulfills all the requirements in the 
lemma and analogous reasonings hold true in the multi- 
partite case where u is replaced by any M-partite target 
state. Moreover, in the bipartite case Prop. pleads imme- 
diately to a simple sufficient distillability criterion, since 
for Gaussian states it is known that Ed > is equivalent 
to a simple inequality for the CM |15| . 

(ii) Entanglement of formation: Continuity of Ep was 
proven in |l8| for finite and in for infinite dimensional 
systems with energy constraint. Super-additivity of Ep is 
a notorious conjecture, which is proven to be equivalent 
to additivity of Ep and to many other additivity con- 
jectur es 1 161. These are known for various special cases 
fcf.|3. ll6lll7| ) but remain to be proven in general. 

(Hi) Squashed entanglement is indeed strongly super- 
additive [20j and its continuity was proven (for finite di- 
mensions) in Ref. |2lj . 

(iv) Logarithmic negativity: En is additive, but fails 
to be strongly super-additive. In fact, En does not only 
fail the requirements for the proposition — Eq.lj^J) turns 



out to be false in this case. As simple counterexample is 
given by the state \tp) = VI - A 2 |00> + A| 11) with A = ± 
In this case E N (ip) ~ 0.57 whereas E N (ip G ) ~ 0.64 [22fl . 
Hence, despite the fact that E^ can easily be calculated, 
it is not a faithful entanglement measure in the sense that 
a Gaussian approximation of a non-Gaussian state could 
lead to an over-estimated amount of entanglement. 

Finally, it is interesting to note that Gaussian states 
not only give a lower but they also provide an upper 
bound to the entanglement if only the CM is known. In 
fact, it is a simple consequence of the maximum entropy 
property that for a given energy (i.e., tr [r] fixed) the 
entanglement is maximized by a Gaussian state 



SECRET KEY DISTILLATION 
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We will now depart from the discussion of entangle- 
ment and see how Lem.^ can be applied to the dis- 
tillation of a classical secret key from quantum states 
under the assumption of collective attacks. Consider 
the case where two parties share m copies of a quan- 
tum state pab and aim at converting these into rm bits 
of a secret key under local operations and public com- 
munication. Allowing for the worst case scenario, in 
which an eavesdropper is given the entire purifying sys- 
tem of pab = ^e[\^ abe) abe\] , this can be under- 
stood as a mapping ^a^ e — > cr® rm <g) p E , where the 
secret bits a = ±(|00)(00| + |H)(H|) are asymptotically 
uncorrelated with the state p E of the eavesdropper. We 
call the asymptotic supremum over all achievable rates 
r the distillable secret key K coVl (pab) of the state. By 
the same reasoning as for the distillable entanglement, 
j£co\\ (together with its multipartite generalizations with 
a = |(|0 ... 0)(0 ... 0| + 11 ... 1)(1 ... 1|)) has the proper- 
ties of being additive and strongly super-additive. Hence, 
under the assumption of continuity Lcm.n implies the 
following. 

Proposition 2 Consider an M-partite system with an 
arbitrary, finite, number of modes per site. Then for ev- 
ery given finite first and second moments the Gaussian 
state pc provides a lower bound to the distillable secret 
key K col \p) > K coll (p G ). 



CHANNEL CAPACITIES 

Let us finally apply Lem. ^ to the task of transmit- 
ting classical information through a Bosonic Gaussian 
quantum channel T 0, 0] . The latter may describe opti- 
cal fibres, harmonic chains or any other Bosonic system 
for which we can describe the evolution in terms of a 
quadratic Hamiltonian H acting on system plus environ- 
ment 



T(p) = tT em \v(p®\<t>)(4>\ em )vi 



V = e 



■Hi 



(12) 
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The classical capacity C of a quantum channel is the 
asymptotically achievable number of classical bits that 
can be reliably transmitted from a sender to a receiver 
per use of the channel. To make this a reasonable concept 
in the infinite dimensional setting, one imposes an energy 
constraint to the encoding. That is, any allowed set of 
input states pi with respective probabilities Pi is such 
that the average state p = X^-PiA satisfies an energy 
constraint p e JC = {p\ V\ tr [(Q] + P j)p] < «}• Under 
this constraint the capacity C(T,K) of the channel is 



Ci(T,/C) 



C(T, JC) 



sup 

{pi,Pi} 



S(T(p)) -J2piS(T( Pi ) 



lim -CAT 

n — >oo ft 



1C 



(13) 



(14) 



Consider now a fixed state p and define p — V(p <8> 
|0)(0|onv)^^. Then we can write 

C 1 (T,p) = S{T(p))-E F (p) . (15) 

If the notorious additivity conjecture is true |l6j . then 
not only Ep satisfies the requirements of Lem.^but also 
C\ = C, i.e., the supremum over all p G JC in Ea. l|15|) 
gives the capacity of T. By Prop. ^ together with the 
maximum entropy property of Gaussian states we have 
then, however, that C(T,p) < C(T, p G ) as both terms 
in Ea.l)15|l become extremal for the Gaussian state p G , 
which has the same CM as p- Since p G JC iff p G G JC this 
shows: 

Proposition 3 Consider a Bosonic Gaussian channel 
acting on a finite number of modes. Then there is an 
optimal encoding, which achieves the classical capacity 
with a Gaussian p. if the capacity is additive. 

For single mode channels for which the optimal p is a 
symmetric two-mode Gaussian state we then even know 
an optimal ensemble {pi, pi} (which is continuous in this 
case). For such two- mode states it has been shown 
that Ep(p) equals the so-called Gaussian entanglement 
of formation |25|. which in turn implies that the optimal 
ensemble consists of coherent states which are distributed 
in phase space according to an appropriate Gaussian dis- 
tribution. 

There are certainly many other applications of Lem. ^ 
proving extremality of Gaussian states in various con- 
texts. Via a state-channel duality one might also ap- 
ply similar techniques to channels and operations instead 
of states (cf. 6J). In fact, recently Gaussian operations 
turned out to be optimal for certain tasks concerning 
classical teleportation and cloning of coherent states (2(| . 
Moreover, it is straight forward to translate the results 
from the Bosonic to the Fermionic world, i.e, in other 
words to work with the CAR instead of the CCR alge- 
bra. In this case a central limit theorem has been proven 
in and in 7] the case of normal fluctuations on an 
arbitrary observable algebra was studied. 
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